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Abstract. Quantum chromodynamics (QCD) phase diagram is usually plotted as 
temperature (T) versus the chemical potential associated with the conserved baryon 
number (pb)- Two fundamental properties of QCD, related to confinement and 
chiral symmetry, allows for two corresponding phase transitions when T and /jb are 
(Jjj ' varied. Theoretically the phase diagram is explored through non-perturbative QCD 

calculations on lattice. The energy scale for the phase diagram (Aqcd ^ 200 MeV) is 
^ ' such that it can be explored experimentally by colliding nuclei at varying beam energies 

(— I . in the laboratory. In this paper we review some aspects of the QCD phase structure 

as explored through the experimental studies using high energy nuclear collisions. 
Specifically, we discuss three observations related to the formation of a strongly coupled 
^ , plasma of quarks and gluons in the collisions, experimental search for the QCD critical 

' point on the phase diagram and freeze-out properties of the hadronic phase. 
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1. QCD Phase Diagram 

Physical systems undergo phase transitions when external parameters such as the 
temperature (T) or a chemical potential (/x) are changed. A phase diagram provides 
intrinsic knowledge on the structure of the matter under study. In other words, it tells us 
how matter organizes itself under external conditions at a given degrees of freedom. The 
theory of strong interactions, QCD, predicts that nuclear matter at high temperature 
and/or density makes a transition from a state where quarks and gluons are confined 
and chiral symmetry is broken to a state where quarks and gluons are de-confined 
and chiral symmetry is restored [1]. QCD has several conserved quantities: baryon 
number (B), electric charge (Q), and strangeness {S). Each of these is associated 
with a chemical potential. As a result, the QCD phase diagram is four-dimensional. 
Hq and Us are relatively small compared to fi^ (baryonic chemical potential) in high 
energy nuclear collisions [2] . The T and /xb are varied in a typical QCD phase diagram 
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Figure 1. (Color online) Typical phase diagram of QCD. See text for details. Figure 
taken from Ref. [3]. 



as shown in Fig. 1 [3]. At high temperature and density the phase is governed by 
quark and gluon degrees of freedom and is commonly referred to as the Quark Gluon 
Plasma (QGP) [4]. At large densities and low temperatures other interesting phases 
related to neutron star [5] and color superconductivity [6] starts to appear. These 
QCD transitions which occurred in the early universe have the right energy scale to 
be accessible by the experiments. The Fig. 1 shows the parts of the phase diagram 
explored by several accelerator based experimental programs. Experimentally, this is 
done by varying the beam energy. Both, T and fi^, vary as the function of the center- 
of-mass energy (-y/snn) [7]. This strategy is being followed by experimental programs 
at Large Hadron Collider (LHC) at CERN, Relativistic Heavy Ion Collider (RHIC) at 
BNL and Super Proton Synchrotron (SPS) at CERN, and will be followed at Facility for 
Antiproton and Ion Research (FAIR) at GSI and Nuclotron-based Ion Collider fAcility 
(NICA) at JINR. Among these experiments, we discuss in section 2, selected results 
from RHIC which provides evidence for the formation of a QGP [8]. 

Theoretically, finite temperature lattice QCD calculations at /iB = suggest a 
cross-over [9] above a temperature, Tc, of about 170 to 190 MeV from a system with 
hadronic degrees of freedom to a QGP [10]. At large /ie, several QCD based calculations 
find the quark-hadron phase transition to be of first order [11]. Going towards the 
smaller /ie region, the point in the QCD phase plane (T vs. /ie) where the first order 
phase transition ends is the QCD Critical Point (CP) [12]. The focus in the coming 
decade would be on attempts to locate the CP both experimentally and theoretically. 
Current theoretical calculations are highly uncertain about the location of the CP. This 
is primarily because the lattice QCD calculations at finite /xb face numerical challenges. 
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The experimental plan (discussed in section 3) is to vary the ^/snn of heavy-ion collisions 
to scan the phase plane and, at each energy, search for signatures of the CP that might 
survive the evolution of the system. In the last section of the review we discuss the 
thermodynamic properties of the hadronic phase. Finally we end with a summary of 
our current understanding of phase diagram and a brief outlook. 

2. Establishing the Partonic Phase at RHIC 

2.1. Strangeness enhancement and formation of a gluon rich plasma 

Enhancement of strange hadron production in high energy heavy-ion collisions [13] due 
to formation of QGP is one of the four classic signatures in this field. The other three 
being, enhanced direct photon [14] and dilepton [15] production and suppression of J/^ 
production [16] in heavy-ion collisions relative to p+p collisions. In a QGP, thermal s 
and s quarks can be produced by gluon-gluon interactions [17]. These interactions could 
occur very rapidly and the s-quark abundance would equilibrate. During hadronisation, 
the s and s quarks from the plasma coalesce to form mesons. Production by this 
process is not suppressed as per the OZI (Okubo-Zweig-Izuka) rule [18]. This, coupled 
with large abundances of strange quarks in the plasma, may lead to a dramatic increase 
in the production of mesons and other strange hadrons relative to non-QGP p-\-p 
coUisions. 

Such predictions of relative enhancement of strange hadron was challenged by 
an alternate idea of canonical suppression of strangeness in small systems being the 
source of strangeness enhancement in A, S and Vt hadrons in high energy heavy-ion 
collisions [19]. The strangeness conservation laws require the production of an s-quark 
for each s-quark in the strong interactions. The main argument in such canonical models 
is that the energy and space time extensions in smaller systems may not be sufficiently 
large. This leads to a suppression of open-strange hadron production in small collision 
systems. These statistical models fit the data reasonably well [20]. These models 
predict two things: (a) strangeness enhancement in nucleus-nucleus collisions, relative 
to p-\-p collisions, should increase with the valence strange quark content of the hadrons 
and (b) the enhancement is predicted to decrease with increasing beam energy [21]. 
Discriminating between the two scenarios, QGP versus the canonical suppression, using 
the experimental data on fC^, A, S and Vt hadrons has been, to some extent, ambiguous. 
Enhancement of 4>{ss) production (zero net-strangeness, hence not subjected to the 
canonical suppression) in heavy-ion collisions relative to p+p collisions would clearly 
indicate the formation of a gluon rich QGP in these coUisions. This would then rule out 
canonical suppression scenario. 

Figure 2 shows the enhancement of strange hadron production at RHIC [22]. Upper 
panel shows the ratio of strange hadron production normalized to (A^part) in nucleus- 
nucleus collisions relative to corresponding results from p-\-p collisions at 200 GeV. The 
results are plotted as a function of (A^part)- K~ , A and H -|- S are found to exhibit an 
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Figure 2. (Color online) Upper panel: The ratio of the yields of , <j), A and S + S 
normalized to average number of participating nucleons ((A'part)) in nucleus-nucleus 
collisions to corresponding yields in inelastic p+p collisions as a function of (A^part) 
at 200 GeV. Lower panel: Same as above for (p mesons in Cu+Cu collisions at 200 
and 62.4 GeV. The p+p collision data at 200 GeV are from STAR experiment and at 
62.4 GeV from ISR [23]. The error bars shown represent the statistical and systematic 
errors added in quadrature. Figure is taken from Rcf. [22]. 



enhancement (value > 1) that increases with the number of strange valence quarks. 
Furthermore, the observed enhancement in these open-strange hadrons increases with 
collision centrality, reaching a maximum for the most central collisions. However, the 
enhancement of (f) meson production shows a deviation from the ordering in terms of 
the number of strange constituent quarks. More explicitly, this enhancement is larger 
than for K~ and A, at the same time being smaller than in case of S + S. Despite 
being different particle types (meson-baryon) and having different masses, the results 
for K~ and A are very similar in the entire centrality region studied. This rules out a 
baryon-meson effect as being the reason for the deviation of (j) mesons from the number 
of strange quark ordering seen in Fig. 2 (upper panel). The observed deviation is also 
not a mass effect as the enhancement in meson production is larger than that in A 
(which has mass close to that of the 0). Further in heavy-ion collisions, the production 
of mesons is not canonically suppressed due to its ss structure. 
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The observed enhancement of meson production then is a clear indication for the 
formation of a dense partonic medium being responsible for the strangeness enhancement 
in Au+Au collisions at 200 GeV. The observed enhancement in meson production 
being related to medium density is further supported by the energy dependence shown 
in the lower panel of Fig. 2. The meson production relative to p+p collisions is 
larger at higher beam energy, a trend opposite to that predicted in canonical models for 
other strange hadrons. In addition measurements have shown that meson production 
is not from the coalescence of KK and is minimally affected by re-scattering effects 
in the medium [24]. Measurements also indicate that (j) mesons are formed from the 
coalescence of seemingly thermalized strange quarks [25]. All these observations put 
together indicate the formation of a dense partonic medium in heavy-ion collisions where 
strange quark production is enhanced. This in turn suggests that the observed centrality 
dependence of the enhancement for other strange hadrons (seen in Fig. 2) is likely to be 
related to the same reasons as in the case of the 4> meson, that it is due to the formation 
of a dense gluon rich partonic medium in the collisions. These experimental data rule out 
the possibility of canonical suppression being the only source of the observed strangeness 
enhancement at -s/snn— 200 GeV. 

2.2. Jet Quenching and highly opaque medium 

One of the most exciting results to date at RHIC is the discovery of suppression in 
the production of high transverse momentum (pr) mesons in nucleus-nucleus collisions 
compared to corresponding data from the binary coUision scaled p-|-p coUisions [26] . This 
has been interpreted in terms of energy loss of partons in QGP. This phenomena is called 
as the jet quenching in a dense partonic matter [27]. The energy loss by energetic partons 
traversing the dense medium formed in high-energy heavy-ion collisions is predicted to 
be proportional to both the initial gluon density [28] and the lifetime of the dense 
matter [29]. The results on high-p^ suppression are usually presented in terms of the 
nuclear modification factor (Raa), defined as: 



where the overlap integral Tab = ^binary/ (^inelastic with Nunary being the number of 
binary collisions commonly estimated from Glauber model calculation [30]. 

In Fig. 3 we show the RHIC data for the Raa{pt), for various produced mesons [31] 
and direct photons [32] in central Au-|-Au collisions at midrapidity. A large suppression 
in high meson production is observed, and those for tt'^'s being almost fiat at 
Raa — 0.2 up to 20 GeV/c. The figure also shows that the level of suppression for 7r°'s, 
77's and (/)-mesons are very similar, which supports the conclusion that the suppression 
occurs in the partonic phase, not in the hadronic phase. This strong suppression of 
meson production is in contrast to the behavior of direct photons, also shown in the 
figure. The direct photons follow binary scaling (i.e. Raa — 1) or no suppression. This 
is a strong evidence that the suppression is not an initial state effect, but a final state 
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Figure 3. (Color online) Compilation of the nuclear modification factor (Raa) for 
mesons and direct photons as measured in RHIC experiments at midrapidity for central 
Au+Au collisions at ^snn = 200 GeV. Also shown are the i?dAu for charged pions for 
y'SNN = 200 GeV. The lines are results from various model calculations. See text for 
more details. 



effect caused by tlie liigli density medium witli color cliarges created in tlie collision. 
This is further consolidated by a demonstration through a controlled experiment using 
deuteron on Au ion collisions, which gave a RdAu{pT) ~ 1 for tt^ at midrapidity and 
high pt [33] . 

The various curves in the Fig. 3 represents different model calculations. The dashed 
curve shows a theoretical prediction using the GLV parton energy loss model [28] . The 
model assumes an initial parton density dN/dy = 800 — 1100, which corresponds to 
an energy density of approximately 5-15 GeV/fm^. The lower dashed curves are for 
higher gluon density. The precision high px data at RHIC has been used to characterize 
medium density fairly accurately. The conclusion being that the medium formed in 
central Au+Au collisions at RHIC has a high degree of opacity [34]. In addition, 
theoretical studies suggest that for a given initial density, the Raa{pt) values are also 
sensitive to the lifetime (r) of dense matter formed in heavy- ion collisions [29]. The 
sohd curves are predictions from Ref. [29] at ^snn = 200 GeV with r = 10 fm/c (i.e. 
larger than the typical system size of ~ 6-7 fm). The parton energy loss calculations 
discussed above attributes the opacity to plasma induced radiation of gluons, much like 
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ordinary bremsstrahlung of photons by electrons. However, the quantitatively large 
suppression pattern observed at high p^, for both light hadrons and those involving 
heavy quarks [35] , showed that the mechanism of energy loss is far from being a settled 
issue, namely, the relative contribution of radiative and collisional forms. As an example, 
shown in Fig. 3 is a comparison of the data to theoretical results (dot-dash curves) on 
Raa from models that consider only coUisional energy loss [36] . This model gives Raa 
values at high px close to the measured values and similar to corresponding values from 
models having only a radiative mechanism for parton energy loss. 

Leading particle measurements, such as the ones shown in Fig. 3 suffer from a 
number of limitations, (a) Leading hadrons come from a mixture of parent quarks and 
gluons. (b) As a fragmentation product, the energy of a leading hadron is not a perfect 
proxy for the energy of the parent parton as it samples a wide range of partonic energies. 
In future we should look forward to potentially three interesting measurements, (i) The 
7-jet process potentially provides access to the underlying scattered parton's energy. 
Measurements of the distribution of particles from the jet opposite, in azimuth, to the 
tagged photon reveals how much energy was lost, and how it was redistributed, by the 
colored parton as it traversed the medium [37] . (ii) Another method is through the full 
reconstruction of jets in heavy- ion collisions [38] . Beyond producing a far better proxy 
for the energy of the parent parton than a leading hadron, this technique allows one 
to trace the evolution of energy flow in directions both longitudinal and transverse to 
the direction of the parent parton. Both these methods arc under active pursual at 
RHIC. (iii) Another important feature of jet quenching is provided by partonic identity. 
While it is difficult to disentangle light quarks from gluons, especially in a heavy-ion 
environment, charm and bottom can be easily tagged by the existence of a charmed 
or bottom hadron in the final state. Due to their large masses the charm and bottom 
quarks are predominantly produced via hard scattering in the initial stage of the high- 
energy heavy-ion collision. The final state spectra can therefore serve as a sensitive 
tool to probe in-medium rescattering and interactions responsible for thermalization. 
This will also allow to study the non-Abelian feature of QCD that results in the gluons 
losing more energy than quarks in the medium [39]. Plans are in place to measure 
the cross-sections and transverse momentum spectra of hadrons with open and hidden 
heavy flavor at RHIC with new detector upgrades. This will also provide useful data to 
understand the different mechanisms of energy loss: collisional versus radiative. 

2.3. Partonic collectivity and low viscosity 

Elliptic flow, V2, is an observable which is thought to reflect the conditions from the 
early stage of the collisions [40]. In non-central heavy-ion collisions, the initial spatial 
anisotropy of the overlap region of the colliding nuclei is transformed into an anisotropy 
in momentum space through interactions between the particles. As the system expands 
it becomes more spherical, thus the driving force quenches itself. Therefore the elliptic 
flow is sensitive to the collision dynamics in the early stages. It is measured, by 
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Figure 4. (Color online) Top left plot: The elliptic flow V2 for (a) light quark hadrons 
and (b) strange quark hadrons. The data is from the minimum bias Au+Au collisions 
at midrapidity for ^/s^m = 200 GeV [41]. Top right plot: Compilation of the number 
of constituent quark scaled V2 as a function of the scaled transverse kinetic energy [42] . 
Bottom plot: Compilation of rj/ s extracted from various measurements in heavy-ion 
collisions at RHIC. 



calculating (cos(2(0 — where is the azimuthal angle of the produced particles 

and ^ is the azimuthal angle of the impact parameter, and angular brackets denote an 
average over many particles and events. 

Figure 4 top left panel shows the RHIC results on f 2 of light quark (tt and p) and 
strange quark (0 and f2) carrying hadrons in Au+Au collisions at 200 GeV [41]. Three 
very distinct experimental observations can be made, (a) At low (< 2 GeV/c) the 
heavier hadrons have smaller V2- Such a mass ordering is expected in hydrodynamics 
calculations of U2(pt) for identified particles [43]. (b) At the intermediate range of 2-5 
GeV/c it is observed that baryons have higher V2 than mesons. The 0-meson V2 plays a 
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crucial role in establishing this baryon-meson difference. Such a separation of baryons 
and mesons in the intermediate pt range has been also observed in measurements of 
the nuclear modification factor, Rcp [44]. These results are consistent with calculations 
from quark recombination models [45] implying the de-confinement of the system prior to 
hadronisation. (c) Comparison between the V2 results for light quark carrying hadrons 
to those from strange quark carrying hadrons indicates both types of hadrons show 
similar magnitude of V2- The multi-strange hadrons (0 and Q) have relatively low 
hadronic interaction cross sections and freeze-out early, hence are considered as the 
most promising probes of the early stages of the collision. All the above results indicates 
that a substantial amount of collectivity has been developed at the partonic stage of the 
heavy-ion collisions. In fact none of the available hadronic models are able to account for 
the observed magnitude of f 2 at ^snn = 200 GeV. Only models which have additional 
partonic interactions explain the (^2) values [46]. 

When elliptic flow V2 is plotted versus transverse kinetic energy (mx — mo), both 
divided by the number of constituent quarks, the V2 for all identified hadrons as well as 
light nuclei below (ttit — ttiq) ~ 1 GeV/c^ falls on a universal curve [47]. mo is the mass 
of the particle. This scaling behavior as shown in Fig.4 further strengthens the evidence 
for formation of partonic matter during the Au -|- Au collision process at 200 GeV. It 
is very hard to explain this observed pattern in a scenario where only hadronic matter 
exists throughout the interaction, whereas the hypothesis of coalescence of hadrons from 
de-confined quarks offers a ready explanation. Turn-off of the scaling at a given beam 
energy would indicate the hadronic side of the phase boundary. 

The V2 measurements of light quark carrying hadrons, the nuclear modification 
factor and V2 for heavy quark carrying hadrons and the differential px correlations 
for charged hadrons have been used to extract information of a dimensionless ratio, 
shear viscosity to entropy (rj/s), for the medium formed in heavy- ion collisions at 
RHIC [48]. Figure 4 bottom panel shows the compilation of 77/s extracted from various 
measurements [49]. It is observed to lie between that conjectured from quantum theory 
{ri/s ~ l/An) and those for liquid Helium at Tc. Such a low value of r]/s (within a factor 
1-10 of the quantum limit) indicates that the matter formed in heavy-ion collisions at 
RHIC has low viscosity, hence is a strongly coupled system. 

3. QCD critical point and thermalization 

The CP is a landmark point in the QCD phase diagram, observation of which will make 
the QCD phase diagram a reality. A close collaboration between the experiments and 
theory perhaps will lead to its discovery. The first step in this process is to establish 
an observable for CP which can be measured experimentally and can be related to 
QCD calculations. In this context, it is important to recall that for a static, infinite 
medium, the correlation length (^) diverges at the CP. ^ is related to various moments 
of the distributions of conserved quantities such as net-baryons, net-charge, and net- 
strangeness [50]. Typically variances (cr^ = {{AN)^); AN — N — M; M is the mean) of 
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these distributions are related to ^ as ~ [51]. Finite size and time effects in heavy- 
ion colhsions put constraints on the values of ^. A theoretical calculation suggests ^ 
^ 2-3 fm for heavy-ion collisions [52]. It was recently shown that higher moments 
of distributions of conserved quantities, measuring deviations from a Gaussian, have a 
sensitivity to CP fluctuations that is better than that of a^, due to a stronger dependence 
on ^ [53]. The numerators in skewness {S = {{AN)^) /a^) goes as ^^'^ and kurtosis [k 
= [((AA^)'*) /a^] - 3) goes as A crossing of the phase boundary can manifest itself 
by a change of sign of 5" as a function of energy density [53, 54]. 
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and baryon chemical potential values at chemical freeze-out extracted from particle 
yields using a thermal model. 
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The first connections between QCD calculations and experiment has been recently 
made [55]. Lattice calculations and QCD-based models show that moments of net- 
baryon distributions are related to baryon number (AA'b) susceptibihties (xb = 
((aa^b) ) _ y .^^ volume) [56, 58]. Then one can construct ratios such as: 
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which do not contain the volume and therefore provide a direct and convenient 
comparison of experiment and theory. In the above expressions the left hand side of 
each equality can be measured in an experiment whereas the right hand side can be 
calculated by lattice QCD. Close to the CP models predict the AA^b distributions to 
be non-Gaussian and susceptibilities to diverge causing the experimental observables to 
have large values. The experimental values should also be compared to those expected 
from statistics, for example if the p and p distributions are individually Poissonian then 
Ka^ for net-protons is unity. 

Experimentally measuring event-by-event net-baryon number is difficult. However, 
the net-proton multiplicity (A^p — A^p = AA^p) distribution is measurable. Theoretical 
calculations have shown that AA^p fluctuations reflect the singularity of the charge and 
baryon number susceptibility as expected at the CP [59]. Non-CP model calculations 
show that the inclusion of other baryons does not add to the sensitivity of the 
observable [55]. 

Figure 5 shows the energy dependence of Sa, kg'^ and ^ for AA^p, compared to 
lattice QCD [56] and Hadron Resonance Gas (HRG) model which does not include a 
CP [57] . The experimental values plotted are for central Au-|-Au collisions for ^snn = 
19.6, 62.4 and 200 GeV. The lattice calculations, which predict a CP around /ie ~ 300 

MeV, arc carried out using two-flavor QCD with number of lattice sites in imaginary 
time to be 6 and mass of pion around 230 MeV [56]. The ratios of the non-linear 
susceptibilities at finite are obtained using Padc approximant resummations of the 
quark number susceptibility series. The freeze-out parameters as a function of y^SNN 
are from [60] and Tc = 175 MeV. 

Prom comparisons of the experimental data to the HRG model and the lack of non- 
monotonic dependence of ku'^ on y^SNN studied, one concludes that there is no indication 
from the current measurements at RHIC for a CP in the region of the phase plane with 
ji-Q < 200 MeV. Although it must be noted that the errors on the experimental data 
point at 19.6 GeV is quite large due to small event statistics. It is difficiilt to rule out 
the existence of CP for the entire region below 200 MeV. The extent to which these 
results can do that is guided by the theoretical work. In addition, the expectation of the 
extent of the critical region in /ib is thought to be about 100 MeV. The results discussed 
here form the baseline for the future CP search program at RHIC [6 1] . However the fact 
that the data shows excellent agreement with HRG and Lattice QCD, both of which 
assume thermalization, is another non-trivial indication of attainment of thermalization 
(some other measurements are discussed in next section) in heavy-ion collisions. Such 
a conclusion is drawn for the first time using fluctuation measurements. 

With the idea that the rise and then fall of observables sensitive to CP as /ib 
increases should allow us to ascertain the {T ^^-q) coordinates of the critical point, the 
beam energy scan program at RHIC has started. The flrst phase of the experimental 
program at RHIC is expected to be completed in 2010-2011. This phase is expected to 
cover a ^/s^ region of 39 to 5.5 GeV, which corresponds to a /xb range of 112 to 550 
MeV [62]. 
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Figure 6. (Color online) Compilation of invariant yield of produced particles at 
midrapidity versus pt in central Au+Au collisions at ^snn = 200 GeV [2, 25, 31, 32, 
63]. The lines show the possibility of having three regions in the spectra where the 
dominant mechanism of particle productions are different. See text for more details. 



4. Hadronic Phase 

The measured hadron spectra reflect the properties of the bulk matter at kinetic freeze- 
out, after elastic collisions among the hadrons have ceased. More direct information 
on the earlier stages can be deduced from the integrated yields of the different hadron 
species, which change only via inelastic collisions. The point in time at which these 
inelastic collisions cease is referred to as chemical freeze-out, which takes place before 
kinetic freeze-out. RHIC experiments have measured the pt distribution of variety 
of particles over a wide range of p^ at midrapidity. A sample of the invariant yield 
( (27rpT^)X'^PT (^^^Z*^)"^) '^^ produced particles at RHIC for central Au+Au collisions at 
y/suN = 200 GeV [2, 25, 31, 32, 63] is shown in Fig. 6. Like the pt dependence of 
V2 discussed in the section 2.3, here also we can separate the spectra into 3 regions 
based on the dominant mechanism of particle production. The low p^ (< 2 GeV/c) 
is explained by thermal model based calculations [64], intermediate px (2-6 GeV/c) 
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by parton recombination based approaches [45] and high px (> 6 GeV/c) by including 
pQCD based processes or jet production [28, 29]. The only statistical distribution which 
so far seems to sucessfully describe the pt spectra and V2{pt) over a wide momentum 
range is the one based on Tsallis statistics [65]. 

In this section we concentrate on low px part of the spectra for rest of the 
discussions. The transverse momentum distributions of the different particles contain 
two components, one random and one collective. The random component can be 
identified as the one that depends on the temperature of the system at kinetic freeze- 
out (Tfo). The collective component, which arises from the matter density gradient 
from the center to the boundary of the fireball created in high-energy nuclear collisions, 
is generated by collective flow in the transverse direction, and is characterized by its 
velocity Pt- Assuming that the system attains thermal equilibrium, the blast-wave 
(BW) formulation [66] can be used to extract Tfo and {/3t) ■ The transverse flow velocity 
of a particle at a distance r from the center of the emission source, as a function of the 
surface velocity {(3s) of the expanding cylinder, is parameterized as Prir) = /3s{r/R)^, 
where n is found by fitting the data. The transverse momentum spectrum is then 

dN , ^ /prsinhp(r)\ f tht cosh. p(r)\ 

where Iq and Ki are modified Bessel functions and p(r) = tanh"^ Prir)- 

Fits to the identified hadrons pr distributions at midrapidity for Au-|-Au coUisions 
at y/SNN — 200 GeV using Eq. 2 are carried out. The extracted model parameters 
characterizing the random (generally interpreted as a kinetic freeze-out temperature 
Tfo) and collective (radial flow velocity (Pt)) are shown in Fig. 7 in terms of confldence 
level (x^) contours, for various impact parameters of the collision. As the collisions 
become more and more central, the bulk of the system, dominated by the yields of 
tTjKjP have lower kinetic freeze-out temperature and develops stronger collective flow. 
On the other hand, even for the most central collisions, the spectra for multi-strange 
particles and fl appear to reflect a higher freeze-out temperature. 

Within a statistical model in thermodynamical equilibrium, the particle abundance 
in a system of volume V can be given by 

'''''' -Wr^^i^;^lE ^^)^/ ^' 

where Nj, is the abundance of particle species i, gi is the spin degeneracy, Bi and Si 
are the baryon number and strangeness number, respectively, Ei is the particle energy, 
and the integral is taken over all momentum space [2]. The model parameters are the 
chemical freeze-out temperature (Tch), the baryon (/xb) and strangeness (fis) chemical 
potentials, and the ad hoc strangeness suppression factor (75). Measured particle ratios 
are used to constrain the values of T^i and fi^ at chemical freeze-out. 

Figure 7 compares STAR measurements of integrated hadron yield ratios for 
central Au-|-Au collisions to statistical model calculations. The excellent agreement 
between data and model are observed. The ratios which include stable and long-lived 
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Figure 7. (Color online) Left plot: The contours, extracted from thermal + radial 
flow fits (without allowance for resonance feed-down), for produced hadrons Tr,K and 
p and multi-strange hadrons </> and il. On the top of the plot, the numerical labels 
indicate the centrality selection. For tt, K and p, 9 centrality bins (from top 5% to 
70-80%) were used for ^/SNN — 200 GeV Au-j-Au colhsions [2]. The results from p-l-p 
collisions are also shown. For and fl, only the most central results are presented. 
Dashed and solid lines are the l-cr and 2-a contours, respectively. Right plot: Ratios 
of pT-integrated mid-rapidity yields for different hadron species measured in STAR 
for central Au-f-Au collisions at ^/snn = 200 GeV. The horizontal bars represent 
statistical model fits to the measured yield ratios for stable and long-lived hadrons. 
The fit parameters are T^h = 163 ± 4 MeV, /iij = 24 ± 4 MeV, 7^ = 0.99 ± 0.07. 
The variation of 7s with centrality is shown in the inset, including the value (leftmost 
point) from fits to yield ratios measured by STAR for 200 GeV p-l-p collisions. 



hadrons through multi-strange baryons, is consistent with the hght flavors, u, d, and 
s, having reached chemical equilibrium (for central and near-central collisions only) at 
Tch = 163 ± 5 MeV. The deviations of the short-lived resonance yields, such as those 
for A* and K* collected near the right side of Fig. 7, from the statistical model fits, 
presumably result from hadronic re-scattering after the chemical freeze-out and needs 
to be further understood. 

The saturation of the strange sector yields, attained for the first time in near-central 
RHIC collisions, is particularly significant. The saturation is indicated quantitatively 
by the value obtained for the non-equilibrium parameter 7., for the strange sector for 
central collisions. The temperature deduced from the fits is essentially equal to the 
critical value for a QGP-to-hadron-gas transition predicted by Lattice QCD [67], but 
is also close to the Hagedorn limit for a hadron resonance gas, predicted without any 
consideration of quark and gluon degrees of freedom [68]. If thermalization is indeed 
achieved by the bulk matter prior to chemical freeze-out, then the deduced value of Tch 
represents a lower limit on that thermalization temperature. 
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Figure 8. (Color online) Temperature vs. baryon chemical potential (hb) from heavy- 
ion collisions at various y^J/vjv [61]. The /ib values shown are estimated at chemical 
freeze-out. The kinetic and chemical freeze-out parameters, extracted using models 
assuming thermal and chemical equilibrium from midrapidity particle ratio and 
spectra measurements in heavy- ion collisions. The range of critical temperatures (Tc) 
of the cross-over quark-hadron phase transition at /iB = and the QCD critical point 
from two different calculations from lattice QCD are also indicated. Model-based 
estimates of the range of initial temperature (Tinitiai) achieved in heavy-ion collisions 
based in part on direct photon data at top RHIC [69] and SPS [70] energies are also 
shown. The range of /j,b to be scanned in the RHIC beam energy scan program 
corresponding to y^Sivjv = 5.5 to 39 GeV as well as experiments at SPS and CBM 
are indicated as shaded ellipse. The solid point around T ^ and ^b — 938 MeV 
represents nuclear matter in the ground state. 



5. Summary and Outlook 

In summary, the current understanding of the QCD phase diagram is depicted in Fig. 8. 
From the QCD calculations on lattice it is now established theoretically that the quark- 
hadron transition at /ig = is a cross-over. The critical temperature for a quark-hadron 
phase transition lies within a range of 170-190 MeV. Most calculations on lattice also 
indicate the existence of QCD critical point for fiB > 160 MeV. The exact location 
of CP is not yet known unambiguously. Two such predictions computed on lattice are 
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shown in Fig. 8 for a of 176 MeV [12]. New distinct signatures have been predicted by 
QCD model calculations to locate the critical point in the phase diagram. The specific 
suggestion for a CP search is to look for non-monotonic variation in the products of 
the higher moments of net-proton and net-charge distributions, which are related to 
susceptibilities, as a function of ^/sNN (or T,ij,b)- At top RHIC energies, fluctuations in 
net-proton numbers are found to be consistent with expectations from HRG and lattice 
QCD calculations [55]. These measurements estabhshed the baseline for the CP search 
in the heavy-ion collision program, put constraints on the location of CP in the QCD 
phase plane and added to the evidences for thermalization in central Au-|-Au collisions 
at RHIC. 

High energy heavy-ion collision experiments have seen distinct signatures which 
suggest that the relevant degrees of freedom in the initial stages of the collisions at top 
RHIC energies are quark and gluons [8]. Three such signatures related to strangeness 
enhancement, jet quenching and partonic collectivity are discussed in this paper. The 
initial temperatures (Tinitiai) achieved at top RHIC and SPS energies are obtained from 
models [71] that explain the direct photon measurements from the PHENIX experiment 
at RHIC [69] and from the WA98 experiment at SPS [70]. Prom these models, which 
assume that thermalization is achieved in the collisions within a time between 0.1-1.2 
fm/c, the Tinitiai extracted is greater than 300 MeV at RHIC and greater than 200 MeV at 
SPS. Further, the understanding of suppression in high pr hadron production in heavy- 
ion collisions relative to p+p collisions at RHIC requires a medium energy density >> 
1 GeV/fm^ (critical energy density from lattice for a phase transition). This also shows 
that the medium has a high degree of opacity to propagation of color charges. In addition 
the measurement of elliptic flow and the observation of number of constituent quark 
scaling demonstrates that substantial collectivity has been developed in the partonic 
phase. The magnitude of the flow across several hadronic species and a small value of 
viscosity to entropy ratio extracted from the data supports the idea of formation of a 
strongly coupled system in the heavy-ion collisions. This then also supports the notion 
of creating a liquid with low viscosity in high energy nuclear collisions [72]. 

The experiments have also measured the temperature at which the inelastic 
collisions ceases (chemical freeze-out) and elastic collisions ceases (kinetic freeze-out). 
These temperatures (as shown in Fig. 8) are extracted from the measured particle ratios 
and transverse momentum distributions using model calculations which assume the 
system is in chemical and thermal equilibrium. 

New experimental programs at RHIC, SPS, FAIR and NIC A facihties have been 
designed to explore a large part of the QCD phase digram, covering a /ib range of 20- 
600 MeV. Whereas the experimental program at LHC (probing the cross over region 
of ~ MeV of the phase diagram) have started to provide an unique opportunity 
to understand the properties of matter governed by quark-gluon degrees of freedom at 
unprecedented high initial temperatures (higher plasma life time) achieved in the Pb-|-Pb 
coUisions at 2.76-5.5 TeV [73]. Both at LHC and RHIC, one speciflc observable that has 
the potential to provide a further understanding of system formed in heavy-ion collision 
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are the dileptons. Theoretically they are from the virtual photons, and are different from 
real photons in having a mass. The dilepton mass opens up a new dimension and can be 
used to study time evolution of the system in heavy-ion collisions. For example, recently 
it has been discussed that virtual photon (dilepton) interferometry provide access to the 
development of collective flow with time [74] . Studying the px dependence of the elliptic 
flow and nuclear modiflcation factor for dileptons for masses corresponding to various 
hadrons and beyond will help us understand partonic collectivity and medium opacity. 
Comparing the spectral functions of resonances decaying to dileptons and hadrons will 
let us know about the medium effects. While the slope of the dilepton px distributions 
will tell us about development of radial flow and provide direct evidence of thermal 
radiation of partonic origin in high energy nuclear collisions [75] . 
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